A unitarity analysis on the 1=0 d wave tttt 
scattering amplitude 



Jian-jun Wang, 1 Z. Y. Zhou and H. Q. Zheng 

Department of Physics, Peking University, Beijing 100871, P. R. China 

February 2, 2008 

We study 1=0 d wave tttt scattering phase shift using a proper unita- 
| rization approach. It is verified that the /2(1270) resonance corresponds 

to a twin pole structure: one on the second sheet, another one on the 
third sheet. Also we find fair agreement on /2(1270) pole's mass and 
width between our results and the PDG values. Besides, our analysis 
£NJ ■ reveals the existence of a virtual state in this channel. Its pole loca- 

tion is determined by fitting experimental data and it is found to be 
in good agreement with the prediction from chiral perturbation theory. 
' Our analysis demonstrates that partial wave amplitudes calculated in 

chiral perturbation theory are reliable in the small \s\ region, contrary 
to some unitarization models. 
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In a series of recent publications we have developed a new dispersion ap- 
proach incorporating unitarity and chiral symmetry, to study scatterings be- 
tween pseudo-Goldstone bosons. The approach is found to be particularly 
useful in revealing low energy pole structure of the scattering amplitudes. It 
is demonstrated that the a meson is essential to adjust chiral perturbation 
theory (xPT) to experimental data. £Q Combining with crossing symmetry, 
the mass and width of the a meson are carefully determined [2] and are found 
in nice agreement with the Roy equation analysis. [Sj By analyzing the LASS 
data jl] it is also demonstrated that the k resonance must exist in I=0,J=l/2 
channel of nK scatterings if prediction on the scattering length in the 
corresponding channel is acceptable. [5] It is also found that, interestingly, in 
1=2, J=0 channel of tttt scattering, a virtual state pole is needed in order to 
fit the experimental data and to correctly reproduce the chiral prediction on 
the scattering length. j2] The existence of such a virtual state pole is actually 
a prediction of current algebra and %PT jHj, hence demonstrating that 
provides a consistent description to the scattering T matrix in the small \s\ 
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region. In this note we will further investigate low energy 1=0 d-wave nir scat- 
terings, using the phase shift data from Ref. ^2]- The results provide further 
evidences to support predictions in the small \s\ region. 

We start from a proper parametrization form, for the partial wave 5* ma- 
trix [2J EJ: 

S phy = Y[S Pi -S cut , (1) 
i 

where S Pi are the simplest S matrices characterizing the isolated singularities 
of S phy on the second sheet in the absence of bound states, which are: 



1. For a virtual state pole located at s = s (0 < s < 4m 
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S{s) = - r ~T\~ > 2 

1 — ip(s)a 



and a is the scattering length: 



(3) 
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The kinematic factor pis) = 

2. For a resonance located at zq (and Zq) on the second sheet, the S matrix 
can be written as, 

R( , = M 2 [z ] -s + ip{s)sG[zo] 
1 1 M 2 [z ] - s - ip(s)sG[z ] ' 1 ' 

where 

M 2 [z ] = Re[z ] + Im[2 ] 



Imy z (z - 4ml)} 



Re[^z (z -4ml)] 

G[z o] = . ^ (5) 

Re[ Jz (z - Ami)} 



According to Eq. ((TJ), each pole's contribution to the scattering phase shift, 
denoted as S Pi (s), is additive and easily calculable using Eqs. (0) and In 
Eq. (Q) S cut can be parameterized in the following simple form, 

S cut = e 2ipf(s) 5 (g) 

ft \ - f ( \ , f ( \ s f lm hf(s'), , 8 f lm R f(s') , 

f(s) = h{s) + f R {s) = - / — — ds + - / — — ds , (7) 

ttjl(s — s)s ttJr(s —s)s 

where L = (— oo, 0] denotes the left hand cut on the negative real axis, R 
denotes cuts from the first inelastic threshold to oo, and the dispersion integral 
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is free from subtraction constants. The discontinuity of f(s) on left or right 
hand cut is expressed as 

lm L , R f(s)^-^--log\S^( S )\ . (8) 

The above parametrization is, however, not directly applicable to the study of 
d-wave scatterings since the following threshold constraint has not yet been 
taken into account: 



ReTu(s) = q^la 1 ; + b^q 2 + 0{q 4 )} , (q = -y/s -4ml), (9) 

To remedy this we first recast Eq. J7J) into a thrice subtracted form with a 
once-subtraction at s = 0, and a twice-subtraction at physical threshold 4m 2 : 



s(s-4m 2 w ) 2 r Imf(s') 

7T Jl+r (s'-s)s'(s'-Aml) 2 ' { } 

where the two subtraction constants, f(4m 2 ) and f'(&m%) are not free, they 
are correlated to pole parameters as dictated by Eq. (JHJ): 



i Pi i 



where 
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f(4ml) = Ylf P 

i 

r^o = £4, (ii) 
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G (16 G 2 m/ + 12 m 2 M 2 — 3 
3 (4m T 2 — M 2 f 



for a resonance, and 



a 3 



f v = -a, K = -— 2 (13) 
12 m n z 

for a virtual state. Since the threshold constraint must be exactly obeyed, we 
in the following will make use of Eq. to replace f(4m 2 ) and /'(4m^) in 
Eq. (jlOj) by corresponding pole parameters. 

According to Eqs. (0), © and (J7J), the phase shift from partial wave S 
matrices for elastic irn scatterings consist of contributions from second-sheet 
poles, left-hand cut integral and right-hand cut integral: 

*(*) = £ W + p(*)/(«) (i4) 
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The right-hand cut represents effects from inelastic thresholds and poles from 
higher sheets. We have: [S] 



lm R f(s) 



4p °l - 

^log 1 - Ap(J2 Pn\T ln (s)\ 2 + ■ ■ ■) 



1 



log 1 -4pIm R T + 4p 2 |T(s)| 2 




(15) 



As an approximation, when only poles are considered, we have the following 
parametrization for inelastic amplitudes: jH] 



with T rn the partial width into channel n, T r the total width of resonance r, 
and p n (s) represents the kinematic factor of n-th channel. Especially in above 
(and hereafter) subscript 1 denotes the 2tt channel, and here R starts from 
47r threshold. It should be stressed that the Eq. (J 16)) only works for narrow 
resonances which is only of limited usage, and furthermore, the narrow width 
approximation only works when the energy region is far from any physical 
threshold. 

For the parametrization scheme described above, it should be realized that 
isolated singularities in a couple channel S matrix can play rather different roles 
depending on which sheet they locate. It is therefore necessary to elaborate 
more about the relation between a physical resonance and its pole structure 
in the couple channel situation. Taking the well established ^(1270) 1=0 2 ++ 
resonance ^Uj for example, the PDG value of its mass and width are M = 
1275.4 ± l.lMeV and T = 185.2t|^MeV. The particle decays predominantly 
into 7T7T with a branching ratio Br("7T7r) = 84.8%. The second largest branching 
ratio is the Aw mode, Br(47r) = 10.2% with sizable error bars. The KK 
mode is less important, with Bt(KK) = 4.6%. Most of the experimental 
results listed in Ref. [TU] are from production experiments using a standard 
Breit-Wigner parametrization form for the <i-wave, which means the detected 
^(1270) resonance is on the sheet smoothly connected to the upper edge of 
the unitarity cut (i.e., the place where experimentalists perform experiments). 
In a more formal language, under for example the two channel approximation, 
the experimentally established resonance is a 3rd sheet pole. However, in the 
present situation the picture of a couple channel resonance is more complicated. 
Indeed, phase shift data in the I,J=0,2 channel exhibit a typical resonance 
structure with a mass around 1270MeV, as shown in fig. |5| Nevertheless it is 
easy to understand from Eq. (0) that the jump of the scattering phase shift at 
around 1270MeV is provided by a second sheet pole. In the present scheme, 
the expression for the phase shift of the elastic scattering amplitude measured 
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Figure 1: An illustrative example for a single resonance's contribution (hence 
without interference between different resonances) to the phase shift and the 
inelasticity parameter above the in-elasticity threshold, drawn using Eqs. (|17|) 
and (PI). 



by experiments in the in-elastic region, according to Eq. fl4|) . reads as 

*M = p>PiWs) + ^P/ R {^|^', (olft-4) (17) 

where V stands for principal value. The inelasticity parameter obeys the fol- 
lowing expression, 

V (s)=exp[-2 Pl (s)lm R f(s)] . (18) 

From Eqs. (J15)) and we know that physically observed narrow states in 
production experiments contribute to the inelasticity parameter. Meanwhile 
it can be verified that their contributions to the phase shift of the elastic 
amplitude are small kinks, rather than rapid jumps of approximately 180° (see 
fig-Efor illustration). In other words, if we could have accurate data on phase 
shift and especially inelasticity, we would be able to get more information on 
poles on different sheets. 2 Nevertheless in practice the data are usually not 
good enough to achieve such a goal. 

The left hand integral appeared in Eq. (J10J) converges very rapidly if we 
approximate S phy in the integrand by S xPT . It has been shown that such an 
approximation provides a satisfactory description to the data at qualitative 
level in exotic channels, and also provides at least a self-consistent description 
in non-exotic channels. [21 C] Therefore we in this paper will also use 
result to evaluate the left hand integral. Owing to the thrice-subtraction, 
the cutoff parameter used when evaluating the left hand integral has only a 
tiny influence, and increasing the cutoff parameter will slightly improve the fit 
quality. Thus the cutoff parameter is practically fixed at oo and and is not 
treated ClS db free parameter. 

2 More precisely speaking, under present scheme we can at best determine the partial 
decay width I\ and the total decay width T tot but no more. 
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The integral on the r.h.s. of Eq. (|10|) behaves as 0(s 2 ) when s is large. 
The coefficient of the 0(s 2 ) term has a negative sign and it has a quite large 
effect at large s (for example, when s > 1.4GeV). Thus it is demanded that the 
polynomial on the r.h.s. of Eq. (|l(Jj) be large enough to counteract the large 
negative contribution from the integral when s is large. It can be shown that 
the second sheet pole / 2 (1270), at around M 2 ~ 1.63 GeV 2 and G ~ 0.13, 
provides a negative contribution to the coefficient of s 2 term. That means 
contributions from other poles are necessary. A virtual state or a resonance 
located in some special region can offer a positive term to the coefficient. For 
the latter case it can be verified numerically that, if denoting the 2nd sheet 
resonance pole position as z Q on the complex s plane, then Re[zo] < 4m 2 . On 
the other hand, a bound state will offer a negative term, and a third sheet 
pole contributes through the right hand cut integral and it also provides a 
negative coefficient. Therefore at least one virtual state or one peculiar second 
sheet resonance should exist in the 1=0 c?-wave amplitude. However numerical 
analysis reveals that the latter case is not stable, the additional resonance will 
collapse and convert to two virtual states. Furthermore, one of the virtual 
state moves towards s = (and hence has a vanishing effect) and another 
remains very close to the location of the virtual pole when only the latter is 
used in the fit. Therefore the conclusion is that a virtual state pole is needed in 
order to correctly describe the data. This observation is remarkable, since the 
virtual pole's contribution to the phase shift is very small and hence is hard 
to be directly seen from the data. On the other hand the rf-wave two-loop 
%PT amplitude does predict a virtual state pole at about 1.95 MeV, and the 
one loop amplitude at 1.52 MeV using the low energy constants provided by 
Ref. [TI]. 

Besides the ^2(1270) resonance as described earlier, Ref. [TU] also provides 
several other 2 ++ resonances below 1.7GeV: ^2(1430), /2(1525), /2(1565), ex- 
cept for the ^2(1270) state. In the present approach they are all contained in 
the right hand integral and contribute to rj as small dips, and to the phase 
shift as small kinks. Indeed the CERN-Munich data^2] as shown in fig. |2] 
may reveal some structure in the region ^Js > 1.4GeV, however the data are 
not accurate enough to be useful in distinguishing those higher mass poles. 
The data on rj only exhibits ^2(1270) clearly, so that fit to 5 and 77 fails to 
reproduce other poles unambiguously. Nevertheless these higher mass poles 
do play a role in reducing the fit total width of fl 11 , meanwhile they play very 
little role in decreasing the total x 2 and altering any other output. In the 
following we perform a 12 parameter fit to the CERN-Munich data on 1=0 d 
wave 7T7T scattering phase shift up to 1.6GeV. The 12 parameters consist of 6 
relevant parameters and 6 'irrelevant' parameters. The former are: 1 for the 
virtual pole position, 2 for /' 2 n (1270) and 3 for i / 2 [n (1270) (mass, total width 
and partial width into tttt). The irrelevant parameters are the total widths and 
BR^irir) of the 3 higher mass states. The extra resonances' masses are fixed 
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Table 1: Phase shift value for the 1=0, J=2 txtx scattering. Error bars are 
obtained using an error matrix generated by 12 fit parameters. 

at their values as provided by Ref. ^U]. The fit results are the following, 

Xlo.f. = 96.5/(80 - 12) ; 
M v = 1.73 ± O.OAMeV ; 

Mf 2 = 1270.6 ± 3.2MeV , = 156.5 ± 5.6MeV ; 

Mf* = 1278.6 ± 22.2MeV , V^(tat) = 225.0 ± 76.2MeV ; 

BR(txix) = 92.1 ±2.4% ; 

a° 2 = (1.89 ± 0.09) x 10~ 3 m; 4 , b° 2 = (-2.00 ± 0.02) x 10~ 4 m; 6 ; 

(19) 

where a® and b® are threshold parameters defined in Eq. (JHj). We also provide 
in tabled a few results of phase shift in order to compare with the results from 
Ref. [T3]. The results on additional resonances are not included in Eq. (fT9l . 
since they totally disagree with their PDG values. 3 Nevertheless it should be 
stressed that outputs as listed above are affected rather little with or without 
these additional resonances included. Especially a® is always in nice agreement 
with its xPT value, and b® is not. The reason why the b® parameter does not 
agree well with the ;\RT prediction may be that the b\ parameter is related 
to higher order of momentum expansions and is therefore more sensitive to 
distant singularities. A numerical test is made and it is found that a distant 
second sheet pole with a very large width can remedy the discrepancy on b® 
with very small influence to other outputs in Eq. (fT9"j). but the situation does 
not correspond to a local x 2 minimum and no meaningful conclusion can be 
drawn whether there is any physics behind this discrepancy. 

The total x 2 given by Eq. (JTHJ) is not very good. From fig. |2] one sees that 
the data fluctuation between 0.8GeV and l.lGeV contributes a lot to the total 
X 2 - It is hard to solely ascribe the discrepancy to the data itself. Because in 
the energy region KK threshold opens and should contribute at least part 
of the discrepancy. The Eq. (fT9"|) is only a narrow width approximation, and 
we know that one condition for the narrow width approximation to work is 

3 Thc main reason for the failure is because a 3rd sheet narrow resonance's contribution 
to the phase shift and inelasticity as shown by fig. ^ are rather tiny and concentrated in a 
small region of ^/s which makes the fit heavily disturbed by the fluctuation of data. In order 
to get correct information about these 3rd sheet poles, it requires very dense data with very 
small error bars. 
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that the energy region is far away from any threshold. Therefore, the use of 
Eq. (|T§j) may also contribute to some of the discrepancy. 

To conclude, within an appropriate parametrization form for the partial 
wave scattering S matrix incorporated with chiral symmetry, we have clearly 
demonstrated that the /2(1270) resonance exhibits a twin pole structure, i.e., 
it generates simultaneously both a 3rd sheet pole and a 2nd sheet pole. The 
two poles locate at similar position, though on different sheets. This is a pic- 
ture exhibited by a typical couple-channel Breit-Wigner resonance: when the 
resonance's coupling to the elastic channel is dominant, two poles on differ- 
ent sheets coexist. In production experiments one sees the '3rd' sheet pole, 
whereas in phase shift analysis the nearly 180° jump is generated by the sec- 
ond sheet pole. The pole parameters of /2(1270) resonance are determined, 
the mass and width of the 3rd sheet pole, though contain rather large error 
bars, are found in agreement with the PDG value. 

It is also worth noticing that our approach, for the first time, uncovers the 
existence of the virtual state from 1=0 d-w&ve ttw scattering data. Actually 
the existence of such a virtual pole may be proved rigorously without rescuing 
to perturbation theory. 4 Our numerical analysis indicates that the location of 
the virtual pole is in good agreement with the %PT prediction. This is very 
remarkable, since the virtual state contributes tiny to the phase shift and it 
can hardly be seen directly from the data. The virtual pole definitely emerges, 
in the numerical analysis, as long as left hand cut integral provides a negative 
contribution when s is large, hence the conclusion is robust. Our work, to the 
best of our knowledge, is the first one in support of the theoretical prediction 
by analyzing experimental data. 

Combining with our previous analysis on the virtual state in 1=2 s-wave 
7T7T scattering [2], we conclude that the perturbative chiral expansion in the 
small \s\ region is trustworthy, at least it provides a self-consistent and sat- 
isfactory description to the data, after embedding into a correctly unitarized 
parametrization form. This is not a trivial statement. For example the Pade 
approximation produces a singularity structure in total disagreement with x?T 
predictions, just not far below the two ir threshold. ^3] Therefore it is mean- 
ingful to reexamine more carefully the structure of the scattering amplitude in 
the small \s\ region. 
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